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Abstract 

The predicted effects of final state interactions such as colour reconnection are investigated by measur- 
ing properties of hadronic decays of W bosons, recorded at a centre-of-mass energy of y/s ~ 182.7 GeV 
in the OPAL detector at LEP. Dependence on the modelling of hadronic W decays is avoided by com- 
paring W + W~ — > qq'qq' events with the non-leptonic component of W + W~ — > qq'tUi events. The 
scaled momentum distribution, its mean value, (x p ), and that of the charged particle multiplicity, 
(n c h), are measured and found to be consistent in the two channels. The measured differences are: 

A(n ch ) = (n^)-2(n^) = +0.7 ± 0.8 ± 0.6 

A(x p ) = (xp) - (a$* v ) = (-0.09 ± 0.09 ± 0.05) x 10~ 2 . 

In addition, measurements of rapidity and thrust are performed for W + W~ — > qq'qq' events. The data 
are described well by standard QCD models and disfavour one model of colour reconnection within 
the Ariadne program. The current implementation of the Ellis-Geiger model of colour reconnection 
is excluded. At the current level of statistical precision no evidence for colour reconnection effects was 
found in the observables studied. The predicted effect of colour reconnection on OPAL measurements 
of M\y is also quantified in the context of models studied. 
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1 Introduction 



Hadronic data in e + e~ collisions can be characterised by event shape distributions and inclusive 
observables such as charged particle multiplicities and momentum spectra. In addition to tests of 
Monte Carlo models, measurement of the properties of the hadronic sector of W + W~ decays allows 
the question of "colour reconnection" [EJ to be addressed experimentally. The decay products of the 
two W decays may have a significant space-time overlap as the separation of their decay vertices 
at LEP2 energies is small compared to characteristic hadronic distance scales. In the fully hadronic 
channel this may lead to new types of final state interactions. Colour reconnection is the general name 
applied to the case where such final state interactions lead to colour flow between the decay products 
of the two W bosons. At present there is general consensus that observable effects of such interactions 
during the perturbative phase are expected to be small 0. In contrast, significant interference in 
the hadronisation process is considered to be a real possibility. With the current knowledge of non- 
perturbative QCD, such interference can be estimated only in the context of specific models Jl] ||]. 
One should be aware that other final state effects such as Bose-Einstein correlations between identical 
bosons may also influence the observed event properties. 



It has been suggested |2|,|4| that simple observable quantities, such as the charged multiplicity in 
restricted rapidity intervals, may be particularly sensitive to the effects of colour reconnection. Later 
studies (l^,|lT| showed that the initial estimates of [|j were incorrect for a variety of reasons []. More 
recently, studies using the Ellis-Geiger model [||,[| suggested that the effect on the inclusive charged 
multiplicity itself may be larger than previously considered and that the mean hadronic multiplicity in 
W + W~ — ► qq'qq' events, (n^), may be as much as 10% smaller than twice the hadronic multiplicity in 
W + W _ — ► c^lVi events, (n:j? ") [12|. The visible effects of such phenomena are expected to manifest 



themselves most clearly for low momentum particles, e.g. as illustrated in [11|. Therefore, studies 
of the fragmentation function, i.e. the distribution of the scaled momentum, x p = p/Eb ea _m, are also 
relevant. The shape of a charged particle multiplicity distribution may be quantified by its dispersion 
(r.m.s.), D, and so D 4q and D qq<v are measured. 

Some earlier estimates of the sensitivity to colour reconnection have been made within the context 
of given models, comparing "reconnection" to "no reconnection" scenarios for W + W~ — * qq'qq' events. 
In general, both the size and sign of any changes are dependent upon the model considered. At the 
expense of a reduction in statistical sensitivity, the dependence on the modelling of single hadronic W 
decays can be avoided by comparing directly the properties of the hadronic part of W + W~ — > qq'tPe 
events with W + W~ — * qq'qq' events. In the current study, the inclusive charged particle multiplicity 
and the fragmentation function are measured for W + W~ — > qq'qq' and the non-leptonic component of 
W + W~ — > qq'ive events. Charged particles associated with the leptonically decaying W are excluded 
from these measurements. The quantities A(n c h) = (?\h) — 2( n ch^)' AD = D 4q — \/2D qq ^ and 
A(xp) = (xp q ) — (xp q ^) are examined. The mean values of the thrust distribution, (1 — T 4q ), and 
the rapidity distribution relative to the thrust axis, (|y 4q |), are measured to characterise the global 
properties of W + W~ — ► qq'qq' events themselves. These are not expected to be particularly sensitive 
to the effects of colour reconnection ||. 

Models of colour reconnection as implemented in the event generators Pythia (Sjostrand-Khoze 
model H), Ariadne (model ||) and Herwig (model @,|l^]) have been used in previous studies to 
assess the sensitivities of the quantities above [O]. The models Ariadne and a "colour octet" variant^] 
of Herwig [jl4|] predicted shifts in {n^) and (x p q ) similar in size to the expected statistical uncertainty 

1 The dominant effect was to have neglected the polarisation of the W. 

2 Merging of partons to form clusters was performed on a nearest neighbour basis, as a partial emulation of the model 
of Reference [g| . 
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on these quantities using the data studied in this letter. However, in the studies of [11] no additional 
retuning of the models was performed after reconnection effects were included. 



In this letter, emphasis is placed on predictions from reconnection models which have been suitably 
tuned to describe Z° data. The predictions from the reconnection models of Ariadne, Sjostrand-Khoze, 
and Ellis-Geiger are considered, where the tuning of these models is summarised in Section |j. The 
version of Herwig including reconnection effects has not yet been tuned by OPAL and is therefore 
not discussed further herein. 



2 Data Selection 



This analysis uses data recorded during 1997 by the OPAL detector, which is described fully else- 
where [15|. The measurement of luminosity is identical to that in [j[6|| . The integrated luminosity used 



in this analysis is 57.21 ± 0.15(stat.) ± 0.20(syst.) pb 1 at a luminosity weighted mean centre-of-mass 
energy of ^fs = 182.68 ± 0.05 GeV M\ 



Events are selected in the W + W — » qq'^ and W + W — > qq'qq' channels using the likelihood 
selections described in ]18|,|19[]. In total, 433 W + W~ — * qq'qq' and 361 W + W~ — > qq'iVe candidates 



are selected. The Monte Carlo models without colour reconnection and the detector simulation are 
identical to those in |19|| . The models of colour reconnection studied are described in Section |j. 

Charged particles may have up to 159 hits in the jet chamber. Tracks used in the analysis are 
required to have a minimum of 40 hits in the | cos 9\ region^ in which at least 80 are possible. At larger 
| cos 9\, the number of hits is required to be greater than 50% of the expected number and also greater 
than 20. Tracks must have a momentum component in the plane perpendicular to the beam axis of 
greater than 0.15 GeV/c, and a measured momentum of less than 100 GeV/c. The extrapolated point 
of closest approach of each track to the interaction point is required to be less than 2 cm in the r-(j) 
plane and less than 25 cm in z. Clusters of energy in the electromagnetic calorimeter were required 
to have a measured energy greater than 0.10 GeV (0.25 GeV) if they occurred in the barrel (endcap) 
region of the detector. 



3 Data Analysis and Correction Procedure 

The analysis of the event properties follows the unfolding procedure described in [plf| . The distributions 
of 

n chi x p, 1 — T and y are corrected for background contamination using a bin- by-bin subtraction 
of the expected background, based on Monte Carlo estimates. Corrections (described below) are 
then applied for finite acceptance and the effects of detector resolution, after which mean values of 
the distributions are calculated. Each observable is evaluated using two samples of e + e _ — > W + W~ 



events generated using the Koralw [20] Monte Carlo program. The first, which includes initial state 
radiation and a full simulation of the OPAL detector, contains only those events which pass the cuts 
applied to the data (detector level). The second does not include initial state radiation or detector 
effects and allows all particles with lifetimes shorter than 3 x 10~ 10 s to decay (hadron level). Both 
samples are generated at the same yfs. Distributions normalised to the number of events at the 
detector and the hadron level are compared to derive bin-by-bin correction factors which are used to 



3 The OPAL coordinate system is defined such that the origin is at the geometric centre of the jet chamber, z is 
parallel to, and has positive sense along, the e~ beam direction, r is the coordinate normal to z, 9 is the polar angle 
with respect to +z and cj> is the azimuthal angle around z. 
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correct the observed distributions of x p , 1 — T and y to the hadron level. In contrast to the other 
observables which are formed from charged particles alone, thrust is measured using both charged 
particles and clusters of electromagnetic energy unassociated with charged particles. 

A bin-by-bin correction procedure is suitable for the quantities above as the effects of finite res- 
olution and acceptance do not cause significant migration (and therefore correlation) between bins. 
Such a method is not readily applicable to multiplicity distributions, due to the large correlations 
between bins. Instead, a matrix correction is used to correct for detector resolution effects, followed 
by a bin-by-bin correction which accounts for the residual effects due to acceptance cuts and initial 



state radiation, as in previous OPAL multiplicity studies |21|. 



The uncorrected multiplicity distributions for the W + W _ candidate events before background 
subtraction are illustrated in Figure [l], together with the predictions of Monte Carlo events without 
colour reconnection, but including detector simulation. The background prediction is the sum of all 



Standard Model processes, as described by the models used in | 22fl . In general, good agreement is seen 
between the data and predictions from the models. Herwig predicts visibly lower charged particle 
multiplicities than the other models, for example almost two units lower than that of Koralw for 
W + W~ —* qq'qq' events. This is a consequence of the poor modelling of the multiplicity difference 
between quark flavours in the OPAL tuning of Herwig to Z° data, and of the negligible b quark 
content of W decays. 



4 Colour Reconnection Models 



Variants of three models of colour reconnection have been investigated. The Sjostrand-Khoze (SK) 
models are based upon the Lund string picture of colour confinement, in which a string is created that 
spans the decay product partons associated with each W. These strings expand from the respective 
decay vertices and subsequently fragment to hadrons. Before this occurs, at most one reconnection is 
allowed between sections of the two strings. The main scenarios considered are called type I and type 
II in analogy to the two types of superconducting vortices which could correspond to colour strings. In 
the SK I scenario, the two colour flux tubes have a lateral extent comparable to hadronic dimensions. 
The probability for reconnection to occur is proportional to the integrated space-time volume over 
which the two tubes cross, where a (free) strength parameter, p, determines the absolute normalisation. 
In the SK II scenario, the strings have infinitesimally small radii and a unit reconnection probability 
upon their first crossing. A third scenario considered, SK II', is similar to SK II but reconnection 
is only allowed to occur at the first string crossing which would reduce the total string length of 
the system. As described in [^1 , the only tuning necessary for these models is to ensure that the 
Jetset hadronisation model gives a good description of Z° data. Therefore, the same parameters 
were used as for the corresponding sample of non-reconnected e + e~ — ► W + W~ events, which for the 
Sjostrand-Khoze models are generated using Pythia. The fractions of W + W~ — > qq'qq' events in 
which reconnection occurs at yfs = 183 GeV are found to be 37.9% for SK I (usingQ p = 0.9), 22.1% 
for SK II and 19.8% for SK II'. 

The second set of two models is incorporated into the Ariadne Monte Carlo program by its 
author. They may be considered as extensions of the earlier mode]0 by Gustafson and Hakkinen 
as both models were implemented using the Ariadne Monte Carlo program and the same criterion 

4 p — 0.4 results in a similar reconnection probability to the SK II model. The value used here increases the sensitivity 
of the data available by enhancing the effects of reconnection predicted by SK I. 

5 In Ql, at most one reconnection was allowed per event and possible reconnections within a single W were not 
implemented. 
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is employed in the reconnection ansatz to determine whether reconnection is allowed. Perturbative 
QCD favours configurations with minimal string length in hadronic Z° decays [24|. When the partons 
of two W bosons are separating and strings formed between them, it is plausible that configurations 
corresponding to a reduced total string length are favoured. In the reconnection model of Ariadne, 
the string length is defined in terms of the A measure, which may be viewed as the rapidity range 
along the string: A = Y^i hi(m? /rn 2 p ), where rrii is the invariant mass of string segment i and m p 
sets a typical hadronic mass scale. Reconnections are allowed, within constraints of colour algebra 
factors, which lead to a reduction in the total A of the system. The first Ariadne model considered, 
referred to herein as AR 2, restricts reconnections to gluons having energies less than Tw, while the 
second Ariadne model, AR 3, does not impose such a restriction. As gluons emitted with energies 
> Tw ~ 2 GeV are perturbative in nature and have been shown to be radiated incoherently by 
the two initial colour dipoles Q, the model AR 3 is disfavoured on theoretical grounds. Multiple 
reconnections per event are permitted and reconnection may occur within different string segments 
of the same W boson. The tuning of the models to describe Z° data is as given in pq| , with the 
following changes to the "a" parameter, MSTJ(41), which governs the hardness of the fragmentation 
function: non-reconnected Ariadne: a = 0.52, AR 2: a = 0.65, AR 3: a = 0.58. The fractions of 
W + W~ — ► qq'qq' events in which reconnection occurs at \/~s = 183 GeV are found to be 51.9% for 
AR 2 and 63.4% for AR 3. 



The third class of models considered is that of Ellis and Geiger || , as implemented in the Vni Monte 
Carlo ||, version 4.12. This model has been used to predict relatively large changes in observables 
such as the mean charged particle multiplicity [p^] , but has so far not been subjected to significant 
comparison with W + W~ data. This model has three variants, called "colour blind" , "colour singlet" 
and "colourful". In the colour blind scenario, colour degrees of freedom are ignored and cluster 
formation from partons proceeds solely on a nearest neighbour basis in space-time. The colour singlet 
case requires the colour degrees of freedom of two partons to add up to a colour singlet before cluster 
formation may take place. In the colourful variant, partons which are not in a relative colour singlet 
state are also allowed to merge with each other, with the net colour degrees of freedom being balanced 
by the emission of additional, coloured partons. Predictions presented using this model correspond 
to the colour blind and colour singlet cases. Predictions from the colourful case are not shown due to 
technical problems. To generate e + e~ — > W + W~ events, with the four fermions generated by Pythia, 
it was necessary to changef] the default value of several program parameters. The default tuning of 
the model is used, with only minor modifications!]. In evaluating predictions from this model, decays 
of Kg , A and ir° were included for consistency with the definition of stable particles given in Section ||] 
(not the default for Vni). Significant energy imbalances were observed in hadronic W + W~ final states 
when comparing the generated particle spectra to the collision energy of 183 GeV: approximately 3% 
of W + W~ — > qq'qq' events differed from 183 GeV by more than 0.5 GeV, while approximately 3% 
of W + W~ — > qq'tPe events had deviations of more than 20 GeV. These imbalances are unrelated to 
changes in the program parameters. 



Results obtained using the Ellis-Geiger model in Vni are found to be rather different to those 
of ||, 12|. Differences predicted between the charged multiplicities in the W + W~ — > qq'qq' and 
W + W~ — ► qq'^v channels are ~2% of {n c ^), with opposite signs in the colour blind and colour singlet 
models. It is noted that with the current version of Vni, the predicted value of (1 — T 4q ) is found to 
be 0.30. This is considerably larger than that of 0.06 given in ||. A small value of (1 — T 4q ) indicates 
that an event has a two-jet like topology, allowing the hadronic showers from neighbouring W bosons 
to overlap, which in general is expected to enhance interconnection effects. 



6 MSTW(3)=2, MSTV(84)=0, 2 (blind, singlet), MSTV(16)=1, MSTV(150)=2. 

7 To ensure that unstable particles are decayed irrespective of their origin, and that charged particles carry a fraction 
of the event energy consistent with naive expectation assuming that pions dominate the final state, two changes are 
made: MSTV(91)=5 and MSTV(94)=0. 
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The Sjdstrand-Khoze, Ariadne and Ellis-Geiger models, after simulation of the OPAL detector, 
are compared with the data distributions of charged particle multiplicity and thrust in the W + W~ — * 
qq'qq' channel. Figure [|(a) presents the same data as in Figure |l|(a) but compares them with the 
predictions of the Pythia and Sjostrand-Khoze models, while Figure |2|(b) compares these data with 
the Ariadne and Ellis-Geiger (colour blind) models. Figure §(c) shows the thrust distribution with a 
subset of these models. The model predictions are represented by smooth curves which pass through 
the centre of the binned predictions of each model. It can be seen that the multiplicity predicted by the 
Ellis-Geiger model is significantly higher than observed in data or predicted by other models studied. 
The same observations are made when comparing this model to W + W~ — ► qq'fz/^ data. The AR 3 
model predicts a charged multiplicity that is approximately 0.7 units lower than the non-reconnected 
Ariadne after detector simulation. From the 1 — T distribution, it is noted that most models provide a 
good description of data, while the Ellis-Geiger model predicts events which are much more spherical. 
With its current implementation and tuning in the Vni program, the Ellis-Geiger model does not 
describe the data and is therefore not used to assess systematics on M\y (see Section fj) . 



5 Systematic Uncertainties 

A number of systematic uncertainties have been considered in the analysis, as summarised in Tables [l] 
and ||. Each systematic uncertainty is taken as a symmetric error and the total uncertainty is defined 
by adding the individual contributions in quadrature. Correlations are explicitly accounted for in 
estimating systematic uncertainties for A(n c h), A(cc p ) and AD. The dependence of the correction pro- 
cedure on the Monte Carlo model is evaluated by comparing results obtained using Koralw, Pythia 
and Excalibur |2(| as the W + W~ signal samples, each using the same tuning of the Jetset ||] 
hadronisation model. The same treatment is applied to estimate model dependence in the subtraction 
of the small ( < O (1%)) W + W~ background contamination in each channel: W + W~ — ► qq'qq' or 
W + W~ — > £ + Vi£'~T'£i events selected as W + W~ — ► qq'iu^, for example. The W + W~ cross-section 
obtained from Gentle [27 was allowed to vary from its central value of 15.71 pb by ±4.3%, corre- 



sponding to the combined statistical and systematic uncertainty on the measured cross-section [ |19[] . 

Hadronisation effects are estimated in a similar manner, unfolding data using Herwig and Ariadne,^ 
and also each of the variants of the Sjostrand-Khoze and Ariadne colour reconnection models. The 
largest variation in any of these seven unfolding tests, dominated by the AR 3 and Herwig models, 
was taken as the systematic uncertainty associated with hadronisation effects. The Ellis-Geiger model 
was not used for unfolding as it does not describe the observed W + W~ data. 

Uncertainties arising from the selection of charged tracks are estimated by varying the track selec- 
tion cuts and repeating the analysis. The maximum allowed values of the distances of closest approach 
to the interaction region in r-<j) and z are varied from 2 cm to 5 cm and from 25 cm to 50 cm, respec- 
tively, and the minimum number of hits on tracks is varied from 20 to 40. The dependence on charged 
track quality cuts is the sum in quadrature of these three effects. 

The dependence on the modelling of the accepted background is the sum in quadrature of two 
components, accounting separately for uncertainties in normalisation and shape. Normalisation de- 
pendence is estimated by taking the largest effect observed when scaling various backgrounds before 
subtraction from the data. The scale factors applied, taken from the uncertainties estimated in [p~9|] , 
are approximately ±10% for the total background in the qq'qq' channel and ±8% in the qq'-ft^ channel, 
and approximately ±11% in Z°/7 — > qq' background alone in both channels. Two-photon backgrounds 

8 The Ariadne events were generated using a different tune of Jetset. 
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are scaled by a factor of 0.5 or 2, while the Z°/7 — > t + t~ background is neglected. Shape dependence 
of the subtracted distribution is estimated by taking the largest variation found when using alternative 
models for a given source of background, such as Herwig in place of Pythia for the Z°/7 process, 
and Pythia in place of Phojet pj| for part of the two-photon process. In the case of the multiplicity 
distributions, the shape dependence was also tested by shifting the multiplicity by ±1 unit. This latter 
variation accommodates the differences between data and models, and the uncertainty on measured 
hadronic multiplicity in Z°/7 — > qq' at y/s — 182.7 GeV [f29jj. 

The four-fermion background was evaluated by comparing the accepted distributions from two 
samples of events generated using the grc4f [^] model: one contains the full set of interfering four- 
fermion diagrams, while the other is restricted to the CC03 set of W pair production diagrams^]. The 
same procedure was followed using Excalibur events. A further test of the four-fermion contribution 
used independent samples of WeI7 e , Z°Z° and Z°e + e~ final states, generated with the grc4f or Pythia 
Monte Carlo models. The systematic uncertainty from the four-fermion background was taken to be 
the largest deviation seen from the results obtained using grc4f. 

Since most of the Monte Carlo samples used in the study were generated at \fs = 183 GeV while 
the data were taken at various centre-of-mass energies between ~ 182 and 184 GeV, the analysis was 
repeated using W + W~ samples generated at 182 GeV and 184 GeV. The larger effect quantifies the 
dependence of the analysis on the choice of centre-of-mass energy. 

To estimate the sensitivity of the results to the unfolding procedure, the dispersions of the mul- 
tiplicity distributions and the mean values (n?), {n^ v ), (xp^), (xp q ), (1 — T 4q ) and (|y 4q |) were 
each evaluated by applying a single correction factor to each of the uncorrected values, rather than 
using the methods described in Section ||. This correction is the ratio of the Koralw prediction 
without detector simulation or initial state radiation, to the corresponding prediction for the same 
observable when these two effects are included. The change in the corrected value gives an estimate 
of the systematic error due to the unfolding process. 

As a final test, the analysis was repeated using additional event selection criteria, based on the 
probability obtained from kinematic fits. These required energy and momentum conservation and 
constrained the masses of the two W boson candidates to be equal. Consistent results were obtained 
and no additional systematic contributions were assigned. 



6 Results 

The measurements of the mean charged particle multiplicities, fragmentation functions, and their 
associated differences, corrected to the hadron level, are: 

(n J) = 39.4 ± 0.5 ± 0.6 

{ri$f v } = 19.3 ± 0.3 ± 0.3 

A(n ch ) = +0.7 ±0.8 ±0.6 

(x 4q ) = (3.16 ± 0.05 ± 0.03) x 10~ 2 

(x^) = (3.25 ± 0.07 ± 0.04) x 10~ 2 

A(x p ) = (-0.09 ± 0.09 ± 0.05) x 10" 2 

9 The W pair production diagrams, i. e. i-channel u a exchange and s-channel Z°/7 exchange, are referred to as "CC03" , 
following the notation of Ref. [Jl3j, p. 11. 
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where in each case the first uncertainty is statistical and the second systematic. Similarly corrected 
measurements of (|y 4q |), (1 — T 4q ) and the dispersions of the charged particle multiplicity distributions, 
are also presented in the upper part of Table ||. 

These results may be compared with the predictions from various Monte Carlo models, with and 
without reconnection effects, given in the lower part of Table [| The predictions of the Ellis-Geiger 
model in Vni are markedly different from data and other models. The mean thrust indicates that the 
events are more spherical than data, the momentum spectrum is significantly softer than data, and 
the mean charged particle multiplicities are significantly higher than in data. Although Vni predicts a 
much softer momentum spectrum and higher charged multiplicity than data, the behaviour is similar 
in both qq'qq' and qq'^z^ channels. 

Good agreement is found between the data and the standard QCD models, with the exception 
of Herwig as discussed in Section ||. The higher value of (?!_?) measured in data than predicted by 
models can be inferred from the uncorrected multiplicity distribution of Figure |l](a), where a shift of 
approximately one unit relative to the corresponding predictions from Koralw and Ariadne may 
be seen. The W + W~ — > qqfi&e data are in good agreement with these models. This difference of 
approximately one unit in multiplicity is reflected in the final A(n c h) obtained for the data. 

It is interesting to note that at the current precision, the mean charged particle multiplicity of 
a single hadronic W decay ((n qq v ) or (n^)/2) is consistent with that of approximately 19.5 units 
predicted for Z°/7 — > qq' at y/s ~ Mw, despite the different flavour composition expected in the two 
cases0. The predicted multiplicity is obtained from a fit [El] to the NLLA QCD calculation for the 
energy evolution of the charged particle multiplicity |3l| , using data between 12 GeV and 133 GeV. 

The difference in mean charged multiplicities in hadronic W decays in qq'qq' and qq'ft^ events, 
A (rich), is found to be consistent with zero at the current level of statistical precision. The dispersions 
in these two channels, which should be related by a factor \/2 in the case where two hadronically 
decaying W bosons in a given event are independent of each other, are also found to be consistent with 
each other. Similarly, the measurements of the mean scaled charged particle momenta are consistent 
in the two channels. In contrast, the colour reconnection models predict non-zero values of A(n cn ) 
and A(xp), as shown in Table [3|. The largest effects are seen for the AR 3 model which is the least 
consistent with the measurements of A(n c h) and A(x p ). 

Figures |3|(a) and (b) show the corrected fragmentation functions for the W + W~ — > qq'qq' and 
W + W~ — * qq'ivi channels, together with predictions from the Koralw, Herwig, AR 3 and Vni 
models. Most models are in good agreement within statistical uncertainties in both cases. An alter- 
native measurement of the mean charged multiplicity may be obtained from the integral of the frag- 
mentation function. The values determined in this way are (n 4q ) = 39.5 ±0.9 and (n qqfe ) = 19.3 ±0.5, 
where the errors given are combined statistical and systematic uncertainties. Figure ||(c) shows the 
ratio of the Xp q distribution to twice the Xp qiu distribution for low particle momenta, x p < 0.2. There 
is no significant discrepancy between any of the viable reconnection models and data. Although the 
behaviour of the AR 3 model indicates differences between the two channels in the low momentum 
region where the effects of colour reconnection are expected to be enhanced, there no indication that 
the data show such a tendency. 

The corrected rapidity distribution for W + W~ — > qq'qq' events, \y\ is shown in Figure together 
with the predictions of various Monte Carlo models. To focus attention on the shape of the distribu- 
tions, the integrals of the model predictions have been normalised to the measured (ra 4q ). All of the 

10 The different flavour composition is predicted to iower the charged particie multipiicity of a singie hadronic W decay 
by ~ 0.6 units relative to Z°/7 — » qq'. 
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viable models provide a good description of the data over the entire range measured. 



7 Estimated effect on Mw measurement 



The effect predicted by each model of colour reconnection on the OPAL measurement of M\y in the 
W + W~ — * qq'qq' channel was investigated by analysing Monte Carlo W + W~ — * qq'qq' events, with 
simulation of the OPAL detector, as if they were data. In each selected event, three kinematic fits are 
performed corresponding to the possible jet-jet pairings, and further criteria are used to select up to 



two reconstructed masses per event, m rec , as described in the accompanying paper 22]. 



For each model a possible shift in the measured mass was estimated in three ways. In the first 
method, a distribution of the event-by-event difference between m rec and mt rue was formed, where 
mtrue is defined as the average of the generated W + and W~ masses. Following the Breit-Wigner 



fitting analysis of 1 22], up to two values of (m rec — "mtrue) may enter the distribution, depending upon 
their relative and absolute kinematic fit probabilities. Events admitted to the distribution are required 
to have radiated less than 0.1 GeV of energy into initial state photons. The shift is estimated from 
the truncated mean, using a range of ±3 GeV. The results of this procedure are given in Table || in 
the column "Method la" . A variant of this method, which accounts for possible biases in the mass 
reconstruction in the non-reconnected samples, defines the colour reconnection effect as the difference 
between the means in the corresponding reconnected and non-reconnected samples. The results are 
given in column "Method lb". 

In the second method, the reweighting analysis used to determine the principal result of [^] was 
applied to each sample of simulated events, giving a fitted mass, Mg t . Assuming that this procedure 
does not introduce any significant bias, the difference between Mst and the input Mw value in the 
models is taken to be the change in measured Mw due to colour reconnection. The results of this 
procedure are given in Table § in column "Method 2a" . The column "Method 2b" presents analogous 
results to Method lb. 



It can be seen that the changes predicted for the non-reconnection models are small and consistent 
with zero. In contrast, a significant difference is visible for the Ariadne models, while effects for the 
Sjostrand-Khoze models are less pronounced. The largest change in Mw is seen in the model AR 3 
which also predicts the most significant change in A(n cri ). This model is disfavoured theoretically 
as noted in Section and experimentally by measurements presented in this letter. Therefore, this 
model is not used to estimate possible systematic effects on M\y measurements. 

The individual biases of the admissible models are used to estimate the influence which colour 
reconnection may have on measurements of Mw- Consistent results are obtained in the three methods. 
Although Method lb is the more statistically precise procedure, the overall systematic uncertainty is 
taken from Method 2b, which accounts explicitly for possible biases inherent in the actual procedure 
with which Mw is extracted from data. The largest individual effect is seen from the AR 2 model, 
with a value of 49 MeV. This is therefore assigned as the current estimate of the systematic uncertainty 
associated with colour reconnection effects on OPAL measurements of Mw- 
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8 Conclusions 



The W + W~ event properties presented here are in good agreement with expectations of standard 
QCD models. A 10% reduction in (nj^) as predicted by Ellis and Geiger |0| is not supported by 
data. This difference in the charged multiplicity is not reproduced by the Monte Carlo program Vni 
in which the Ellis-Geiger model is implemented. The current implementation of the Ellis-Geiger model 
is excluded by the measured event shape data and therefore is not used in estimating the systematic 
uncertainty on M\y from colour reconnection. 

Studies of reconnection phenomena implemented with the Ariadne and Pythia models show that 
changes up to approximately 3% may be expected in (nJJ) and (xp-). Defining A(n c h) and A(x p ) 
using data alone provides a model independent test of possible reconnection effects. The maximum 
shifts in these variables predicted by the models considered are at the level 1-2 standard deviations 
for the current data. There is no indication of the effects of colour reconnection on these observables 
at the current level of statistical precision. 

The effect of colour reconnection on measurements of M\y is estimated using the Sjdstrand-Khoze 
and Ariadne Monte Carlo models of colour reconnection which are theoretically admissible and 
consistent with the currently available data. A model dependent estimate of AMw = ±49 MeV is 
currently assigned to the possible influence that colour reconnection may have on il% measurements 
in the W + W~ — * qq'qq' channel. 
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Systematic variation 


\ ch/ 


\ ch / 


A(n rh ) 


D 4q 


J)qq£u 


AD 


W + W~ model 


0.14 


0.03 


0.13 


0.13 


0.08 


0.16 


Hadronisation model 


0.45 


0.21 


0.41 


0.45 


0.31 


0.11 


Track quality cuts 


0.25 


0.12 


0.25 


0.15 


0.11 


0.08 


Background model 


0.23 


0.05 


0.17 


0.03 


0.02 


0.02 


Four-fermion background 


0.18 


0.05 


0.08 


0.07 


0.02 


0.06 


Centre-of-mass energy 


0.01 


0.11 


0.21 


0.06 


0.06 


0.08 


Unfolding procedure 


0.12 


0.03 


0.17 


0.07 


0.06 


0.10 


Total 


0.62 


0.28 


0.60 


0.51 


0.35 


0.25 



Table 1 : Individual contributions to the systematic uncertainties on charged multiplicity related quan- 
tities. 



Systematic variation 


(4 q ) 




A(x p ) 


(1 - T^) 


(|y 4q l> 




xlO 2 


xlO 2 


xlO 2 


xlO 2 


xlO 2 


W+W~ model 


0.012 


0.008 


0.018 


0.11 


0.28 


Hadronisation model 


0.012 


0.015 


0.019 


0.29 


0.81 


Track quality cuts 


0.019 


0.018 


0.015 


0.11 


0.31 


Background model 


0.003 


0.006 


0.005 


0.27 


0.96 


Four-fermion background 


0.011 


0.019 


0.027 


0.08 


0.23 


Centre-of-mass energy 


0.012 


0.019 


0.018 


0.36 


0.38 


Unfolding procedure 


0.002 


0.019 


0.022 


0.68 


0.17 


Total 


0.031 


0.042 


0.049 


0.87 


1.4 



Table 2: Individual contributions to the systematic uncertainties on the average event properties. 
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A{x p ) 
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xlO 2 


xlO 2 


xlO 2 








Data 




39.4 


19.3 


±0.7 


8.8 


6.1 


±0.2 


3.16 


3.25 


-0.09 


0.240 




1.011 


stat. 




±0.5 


±0.3 


±0.8 


±0.3 


±0.3 


±0.5 


±0.05 


±0.07 


±0.09 


±0.015 




±0.014 


syst. 




±0.6 


±0.3 


±0.6 


±0.5 


±0.4 


±0.3 


±0.03 


±0.04 


±0.05 


±0.009 




±0.014 


KORALW 


0.0 


38.78 


19.38 


±0.01 


8.47 


6.00 


-0.01 


3.21 


3.21 


0.00 


0.240 




1.009 


Herwig 


0.0 


37.24 


18.63 


-0.02 


8.65 


6.14 


-0.03 


3.31 


3.30 


0.01 


0.239 




1.009 


Pythia 


0.0 


38.77 


19.39 


-0.01 


8.49 


5.97 


±0.05 


3.21 


3.21 


0.00 


0.240 




1.009 


SK I 


37.9 


38.39 





-0.38 


8.45 





±0.01 


3.25 





±0.04 


0.239 




1.015 


SK II 


22.2 


38.59 





-0.19 


8.45 





±0.01 


3.23 





±0.02 


0.240 




1.010 


SK II' 


19.8 


38.54 





-0.24 


8.41 





-0.03 


3.23 





±0.02 


0.240 




1.011 


Ariadne 


0.0 


38.47 


19.22 


±0.04 


8.23 


5.82 


0.00 


3.22 


3.22 


0.00 


0.240 




1.009 


AR 2 


51.9 


38.31 


19.31 


-0.31 


8.01 


5.67 


-0.01 


3.24 


3.21 


±0.03 


0.239 




1.013 


AR 3 


63.4 


37.45 


19.17 


-0.90 


7.96 


5.67 


-0.05 


3.31 


3.23 


±0.08 


0.238 




1.024 


VNi(blind) 


> 


68.5 


35.0 


-1.6 


14.8 


10.6 


-0.1 


1.80 


1.74 


-0.06 


0.30 




0.81 


VNl(singlet) 


> 


70.3 


34.4 


±1.4 


15.5 


10.7 


±0.5 


1.76 


1.89 


-0.13 


0.31 




0.80 



Table 3: Summary of measurements performed, after correction to the hadron level. The predictions from a variety of Monte Carlo models, with 
and without colour reconnection effects, are also given. The fraction of events in which at least one reconnection has occurred in W + W~ — > qq'qq' 
events is indicated by the column, p rcco . Statistical uncertainties in the Monte Carlo predictions are typically 1-2 units in the least significant 
digit. In the case of the Sjostrand-Khoze models which do not include reconnection for W + W~ — ► qq'tPe events, the predictions from the 
corresponding "no reconnection" sample (Pythia) were used to determine the quantities A(n c h), AD and A(x p ). 







M w bias (MeV) 




Model 


Method la 


Method lb 


Method 2a 


Method 2b 


KORALW 


-11 ± 11 




-12 + 18 




Herwig 


-13 ±9 




-26 + 15 




Pythia 


+2 ±9 




+30 ± 18 




SK I 


+43 ±9 


+41 ± 13 


+41 ± 18 


+11 + 25 


SK II 


+1±9 


-1 ± 13 


+7 + 18 


-23 ± 25 


SK II' 


+4 ±9 


+2 + 13 


+8 + 18 


-22 + 25 


Ariadne 


-10 ±5 




-7 + 10 




AR 2 


+37 ±5 


+47 + 7 


+42 ± 10 


+49 ± 14 


AR 3 


+100 ± 10 


+110 + 11 


+138 + 18 


+145 + 21 



Table 4: Predictions for the bias in the measured M\y from a variety of Monte Carlo models, with 
and without colour reconnection effects. Simulation of the OPAL detector is included. 
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Figure 1: Uncorrected charged multiplicity distributions for (a) W + W~ — ► qq'qq' events and (b) 
the hadronic part of W + W~ — ► qq?£ve events. Points indicate the data, smooth curves show the 
expected sum of signal and background contributions for a variety of signal models, and the hatched 
histogram shows the expected background. Predictions of Koralw, Pythia and Excalibur are 
indistinguishable from one another. 
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Figure 2: Uncorrected charged multiplicity distributions for W + W~ — > qq'qq' events, compared with 
(a) Pythia and the three Sjostrand-Khoze models of colour reconnection, and (b) Ariadne and 
Ellis-Geiger (colour blind) models of colour reconnection. Similarly, (c) compares the uncorrected 
thrust distribution for W + W _ — > qq'qq' events with a selection of the models, with and without 
colour reconnection. Points indicate the data, smooth curves show the expected sum of signal and 
background contributions for a variety of signal models, and the hatched histogram shows the expected 
background. 
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Figure 3: Corrected x p distributions for (a) W + W~ — ► qq'qq' events, (b) the hadronic part of 
W + W~ — > qq'-fi^ events and (c) the ratio of the W + W _ — > qq'qq' distribution to twice the 
W + W~ — > qjtfi&g distribution. Points indicate the data, with statistical (horizontal bars) and sys- 
tematic uncertainties added in quadrature. Point to point correlations exist in the systematic un- 
certainties. Predictions of various Monte Carlo models (with and without colour reconnection) are 
shown as smooth curves. Vni predictions correspond to the colour singlet variant of the Ellis-Geiger 
model. Sjostrand-Khoze and AR 2 models are not shown as they are essentially indistinguishable from 
Koralw. 
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Figure 4: Corrected rapidity distribution for W + W~ — > qq'qq' events. Points indicate the data, with 
statistical (horizontal bars) and systematic uncertainties added in quadrature. Point to point correla- 
tions exist in the systematic uncertainties. Predictions of various Monte Carlo models (with and with- 
out colour reconnection) are shown as smooth curves. Vni predictions correspond to the colour singlet 
variant of the Ellis-Geiger model. All models are normalised to the measured (n c ^). Sjostrand-Khoze 
and AR 2 models are not shown as they are essentially indistinguishable from Koralw. 
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